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Abstract 

A  detailed  chemical  kinetics  model  has  been  built  to  examine  the  gas-phase  chemistry 
between  isocyanic  acid  (HNCO),  white  fuming  nitric  acid  (WFNA),  N20,  C02  and  water. 
This  kinetics  model  is  able  to  explain  the  gas-phase  ignition  observed  during  hypergolic 
ignition  of  the  ionic  liquid;  l-butyl-3-methyl-imidazolium  dicyanamide  with  WFNA. 
Sensitivity  analyses  have  been  performed  to  examine  the  reaction  pathways  for  ignition. 
Ignition  is  predicted  to  occur  via  an  exothermic  reaction  between  isocyanic  acid  (HNCO)  and 
nitric  acid  (H0N02),  and  subsequent  H0N02  thermal  decomposition  that  has  N02  and  OH 
radicals  as  the  primary  chain  carriers.  A  detailed  understanding  of  the  initiation  processes  in 
the  liquid  phase  is  needed  as  the  1 -butyl-3 -methyl-imidazolium  dicyanamide  and  WFNA 
begin  to  react  to  produce  the  above  pre-ignition  species  for  the  proposed  chemical  kinetics 
model  to  describe  the  ignition  behavior  of  the  system. 
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1.  Introduction 


Meeting  energy  demands  while  minimizing  environmental  impact  is  a  major  challenge  for  the 
21st  century.  Pollution  and  greenhouse  effects  mostly  due  to  the  combustion  of  fossil  fuels 
and  disappearance  of  easily  accessible  fossil  fuel  reservoirs  and  other  combustible  materials 
have  led  to  the  search  for  new  fuels  that  promise  to  be  environmentally  benign  as  well  as  have 
superior  energy  content.  Some  of  these  new  fuel  materials  may  not  be  totally  “green”  but 
they  address,  at  least  partially,  some  of  the  drawbacks  of  the  current  fuels  and  combustibles. 
This  ecological  emphasis  is  not  limited  to  auto  fuels  but  extends  to  all  fuels  and  combustibles, 
including  space  propulsion  fuels.  One  of  the  most  utilized  bipropellant  combinations  in  the 
propulsion  of  satellites  and  interplanetary  vehicles,  mono  methylhydrazine  (MMH)/nitrogen 
tetroxide  (NTO)  falls  under  this  scope.  These  propellants  are  earth-storable  liquids  (non- 
cryogenic),  but  their  vapors  are  carcinogenic  and  toxic.  MMH’s  positive  enthalpy  of 
formation  enables  it  to  sustain  a  decomposition  flame,  and  therefore  can  be  considered  to  be 
more  hazardous  than  conventional  fuels  which  have  negative  enthalpies  of  formation. 
Furthermore,  MMH  compatibility  is  often  very  poor  with  many  materials,  and  upon  its 
contact  with  these,  it  can  rapidly  decompose  exothermically.  Similarly,  NTO  is  also  known  to 
be  corrosive  to  many  materials.  It  would  appear  from  the  above  facts  that  there  are  numerous 
valid  practical  reasons  to  replace  MMH  by  other  fuels  that,  at  least,  reduce  the  risks  in 
handling,  storing  and  transporting  the  fuel.  Finding  NTO  replacements  is  also  of  interest  and 
under  certain  circumstances  the  use  of  nitric  acid  or  other  suitable  oxidants  may  be 
advantageous.  Among  some  of  the  alternatives  proposed  to  replace  MMH  are  ionic  liquids 
(ILs).  ILs  exhibit  low  vapor  pressures,  an  advantage  compared  to  MMH  because  this  ensures 
limited  explosion  risks  and  personnel  exposure  to  potentially  toxic  and  carcinogenic  vapors. 
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The  most  interesting  feature  of  MMH  in  combination  with  NTO  is  that  it  is  hypergolic,  that  is, 
MMH  and  NTO  form  mixtures  that  auto-ignite  at  room  temperature  and  at  low  pressures. 
This  characteristic  is  known  to  provide  high  reliability  in  engine  startups  for  numerous 
satellite  propulsion  systems  since  no  external  devices  (such  as  ignition  sparks)  are  required  to 
ignite  the  MMH/NTO  combination.  Over  the  last  15  years,  the  theoretical  understanding  of 
hypergolicity  in  MMH/NTO  has  sufficiently  matured  to  allow  for  reasonable  predictions  in 
the  ignition  delay  time,  a  quantity  of  significant  importance  in  the  design  of  hypergolic 
engines.  [1-51  The  prospect  of  similarly  predicting  ignition  delay  times  for  recently 
discovered  and  new  [6,7]  ionic  liquid  hypergols  will  require  a  detailed  understanding  of  the 
underlying  chemical  pathways  and  the  reaction  mechanisms  involved  when  ILs  are  treated 
with  NTO  or  any  other  suitable  oxidizer  such  as  red  fuming  nitric  acid  (RFNA),  white  fuming 
nitric  acid  (WFNA)  or  even  hydrogen  peroxide.  Comprehensive  gas-phase  thermodynamics 
and  chemical  kinetics  databases  will  need  to  be  developed  since  the  literature  is  sparse  in  this 
area. 

1.1.  General  description  of  ILs 

ILs  are  synthesized  by  the  reaction  of  a  base  with  a  strong  acid.  ILs  are  therefore  salts  and 
consist  of  cations  and  anions.  Each  ion  and  counter  ion  in  the  IL  will  have  its  own  unique 
chemical  and  physical  attributes,  and  in  a  combination,  the  choice  of  cation-anion  pairs  can 
lead  to  approximately  10 18  different  ionic  liquids  having,  in  principle,  any  desirable  or  tunable 
property.  The  number  of  ILs  reported  in  the  literature  has  expanded  almost  exponentially 
over  the  past  two  decades,  as  have  their  applications  ranging  from  use  in  electro-deposition  of 
metals,  fuel  desulphurization,  propellants,  fuel  cells,  etc. [8]  Theoretical  tools  have  been 
recently  developed  to  understand  some  of  the  ILs’  physical  and  thermodynamical  properties 
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to  better  target  the  utility  of  known  ILs  for  the  given  industrial  process  under  consideration. 
Similarly,  new  simulation  tools  will  need  to  be  developed  to  predict  those  physicochemical 
properties  that  allow  for  an  efficient  approach  to  the  screening  and  designing  of  ILs  for 
hypergolic  applications. 

1.2.  Requirements  for  energetic  ILs 

In  propulsion  applications,  ILs  could  have  several  advantages  compared  to  the  existing 
hypergolic  fuel,  hydrazine  (N2H4)  and  its  methylated  derivatives  MMH  and  UDMH 
(unsymmetrical  dime  thy  lhydrazine,  (CLLLNNLL)).  ILs  have  been  regarded  to  be  thermally 
stable,  and  shown  to  have  negligibly  low  vapor  pressures, [9]  high  densities  and  large  liquid 
ranges.  [10]  For  energetic  ILs,  the  triazolium  and  tetrazolium  cations  and  their  corresponding 
substituted  analogues  have  been  preferred.  The  anionic  counterparts  are  generally  the  nitrate 
(NCV),  the  perchlorate  (CIO4)  or  the  dinitramide  (NfNChhd  ions.  For  hypergolic 
applications,  additional  anions  that  exhibit  high  reactivity  with  oxidizers  need  to  be 
considered.  [71  An  exhaustive  bibliographic  search  shows  that  physiochemical  property  data 
exists  only  for  a  relatively  few  compounds,  even  for  such  basic  properties  such  as  densities 
and  melting  points.  The  preferred  performance  requirements  for  energetic  ILs  in  terms  of 
melting  point,  surface  tension,  density,  viscosity,  etc.  have  recently  been  reported  by  Smiglak 
et  al.  [10] 

1.3.  Ignition  and  combustion  of  ILs 

Jones  et  al.  [11]  have  reported  self-sustained  burning  of  1  -ethyl-4, 5-dimethyltetrazolium 
tetranitratoaluminate  by  thermal  heating  of  the  ionic  liquid  to  about  200°C.  An  ignition  delay 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


4 


of  a  few  seconds  was  observed.  As  a  monopropellant,  this  IL  is  reported  to  have  a  better 
performance  measure  than  hydrazine.  [11]  Li  and  Litzinger  [12]  have  observed  laser-driven 
self-ignition  of  lH-4-amino-l,2,4-triazolium  nitrate  vapors  in  air  at  atmospheric  pressure. 
Ignition  occurred  when  the  heat  flux  reached  -150  W  cm'2,  and  in  helium,  ignition  occurred  at 
a  flux  of  ~200  W  cm"2.  Images  of  the  reported  ignition  in  air  show  that  the  solid  ionic  liquid 
first  melts  with  bubbling,  then  vaporizes  (“smoke”  observed)  and  finally  ignites  in  the  gas 
phase.  Alfano  et  al.  [13]  observed  laser  ignition  of  HAN  and  HEHN  propellant  mixture. 
Recently,  papers  have  been  reported  on  the  potential  use  of  ILs  as  hypergols.  [6,7,14-16] 

1.4.  Chemical  nature  of  vapor  above  ILs 

The  chemical  nature  of  the  vapor  above  an  ionic  liquid  is  difficult  to  determine.  It  is, 
however,  important  to  know  how  the  vapor  pressure  forms  from  an  ionic  liquid  because 
ignition  occurs  in  the  gas  phase  above  the  liquid.  Experimental  findings  and  theoretical 
calculations  have  been  reported  in  the  literature  on  this  issue.  [17-23]  Ionic  liquids  can  be 
divided  into  two  classes:  protic  ionic  liquids  and  aprotic  ionic  liquids.  Protic  ILs  have 
hydrogen  atoms  on  the  ring  nitrogens,  whereas  aprotic  ILs  have  alkyl  groups  on  the  ring 
nitrogens  instead  of  hydrogen  atoms.  For  a  protic  ionic  liquid  abbreviated  as  BH+  X",  for 
instance,  there  could  be  ion  pairs  (BH+  X"(g)),  separated  ions  (BH+(g)  +  X"(g)),  neutral  pairs 
BHX(g)  or  separated  neutrals  (B(g)  +  HX(g))  in  equilibrium  above  the  liquid  surface. 
Apparently,  protic  ILs  only  evaporate  as  separated  neutrals  [23]  whereas  many  aprotic  ILs 
evaporate  as  ion  pairs.  [24]  For  monopropellant  systems  where  a  single  ionic  liquid  is  heated 
up  to  ignition,  a  detailed  understanding  of  the  relative  importance  of  these  processes  would  be 
of  significant  importance  to  properly  model  the  behavior.  On  the  other  hand,  for  room 
temperature  ILs/WFNA  mixtures,  the  gas-phase  composition  may  be  more  complex  because 
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exothermic  reactions  occurring  in  the  liquid  phase  will  produce  additional  gaseous  species  not 
observed  above  the  pure  ionic  liquid.  The  concentrations  of  these  additional  species  likely 
exceed  that  of  the  pure  IL  components,  and  may  ultimately  be  responsible  for  the 
hypergolicity  of  the  system. 

2.  Thermochemistry  of  compounds  formed  during  IL-WFNA  reactivity 

The  most  important  feature  of  the  MMH/NTO  system  is  its  hypergolicity,  which  is  its  auto¬ 
ignition  at  ambient  temperatures  and  pressures.  Recently,  hypergolic  ionic  liquids  have  been 
cited  in  the  literature,  [6,7,  14-161  which  include  mechanistic  information  on  a  class  of  ILs 
which  exhibit  hypergolicity  at  room  temperature  with  WFNA.  The  proposed  mechanism  of 
the  reactions  in  the  liquid  phase  was  used  to  explain  the  FTIR  observations  of  the  gas  phase 
species  above  the  reacting  IL  and  WFNA  mixtures.  In  the  l-butyl-3-methylimidazolium 
dicyanamide  reaction  with  WFNA,  gas-phase  HNO3,  CO2,  N20,  HNCO,  FLO,  and  N02  was 
observed  prior  to  ignition.  HNCO,  N20  and  C02  were  presumably  formed  from 
dinitrobiuret,[25]  N02NHC0NHC0NHN02  (abbreviated  as  DNB,  Fig.  1),  which  decomposes 
thermally  to  give,  intermediate  #1,  N02NHC0NC0  (Fig.  1),  and  intermediate  #2,  NH2N02 
(Fig.  1)  according  to: 
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Fig.  1:  Thermal  decomposition  of  DNB  in  the  liquid  phase  leading  to  two  condensed 
phase  products;  intermediates  #1  and  #2.  The  thermal  decomposition  of  these  products 
is  considered  to  be  the  source  of  the  species  observed  in  the  gas  phase  before  ignition. 

The  two  intermediates  in  the  liquid  phase  can  decompose  to  give  HNCO  +  N2O  +  CO2,  and 
N20  +  H2O,  respectively.  These  compounds,  together  with  HNO3  and  NO2  from  WFNA,  are 
thought  to  lead  to  ignition.  DNB  and  intermediates  #1  and  #2  were  not  observed  in  the  gas 
phase.  This  is  not  surprising  because  DNB  and  the  intermediates  are  highly  reactive  and 
therefore  short-lived  and  are  not  likely  to  be  detected.  DNB  is  known  to  decompose  around 
373  K.[25]  Here,  we  will  attempt  to  explain  the  ignition  by  considering  only  those  species 
actually  observed  in  the  gas  phase.  However,  for  completeness,  the  thermochemistries  of 
DNB  and  the  intermediates  are  also  briefly  discussed.  DNB,  intermediates  #1  and  #2  are  nitro 
compounds  and  their  gas-phase  thermochemistries  can  be  calculated  by  using  the  theoretical 
method  previously  established  for  nitro  and  nitrate  compounds.  [26]  Thermodynamics  data 
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(standard  enthalpy  of  formation  at  298  K,  Cp°  against  T  and  S°  against  T)  on  these  species  are 
needed  for  further  refinements  of  the  derived  gas-phase  ignition  models. 

The  gas-phase  thermochemistry  of  DNB  can  be  calculated  by  using  a  semi-empirical  method 
established  for  nitro  compounds.  Validation  of  this  method  was  presented  by  Osmont  et 
al.  [221  This  procedure  leads  to  a  gas-phase  standard  enthalpy  of  formation  of  DNB  at  298  K 
of  -47.9  kcal  mol"1.  For  intermediate  #  2,  the  nitramide  molecule,  the  data  available  are 
somewhat  inconsistent.  Burcat  and  Ruscic  [27]  in  the  latest  version  of  their  thermodynamics 
database  report  two  quite  different  values  for  the  gas-phase  standard  enthalpy  of  formation  at 
298  K:  -6.2±2.4  and  -0.7±X.X  kcal  mol"1.  This  method  applied  to  NFBNCB  leads  to  +1.4±3.0 
kcal  mol"1,  which  is  more  consistent  with  the  second  value  reported  above.  Similarly,  for 
intermediate  #1,  a  value  of  -44.5±3.0  kcal  mol"1  is  obtained.  The  enthalpy  of  the  reaction 
DNB— >#1+  #2  is  +4.8  kcal  mol"1. 


The  main  species  in  the  gas  phase  above  the  reacting  liquids  are  HNCO  (isocyanic  acid), 
nitric  acid  (HNO3  or  HONO2),  NO2,  N20,  FEO  and  CO2.  The  thermochemistries  of  these 
small  species  are  generally  known  and  are  reported  in  the  CHEMKIN  format  by  Burcat  and 
Ruscic.  [27] 

3.  Detailed  chemical  kinetics  model 

3.1.  Approach 
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It  seems  clear,  from  the  experimental  studies  presented  in  Section  2,  that  gaseous  mixtures  of 
HNCO,  HNO3,  N20,  NO2,  CO2  and  H20  are  able  to  sustain  a  flame.  However,  the  chemical 
pathways  involved  in  ignition  need  to  be  elucidated,  as  they  are  critical  to  the  detailed 
construction  of  a  chemical  kinetics  model.  Generally,  detailed  chemical  kinetics  models  are 
validated  with  experimental  data  obtained  with  premixed  mixtures  at  known  temperatures  and 
pressures,  which  include  ignition  delays,  species  profiles,  fundamental  flame  speeds,  etc. 
obtained  under  “ideal”  laboratory  setups.  In  the  case  of  hypergolic  systems,  as  is  the  case 
here,  it  is  much  more  difficult  to  determine  the  exact  composition,  pressure  and  temperature 
of  the  gas  phase  prior  to  ignition  and  therefore  will  be  much  more  difficult  to  validate  a  given 
chemical  kinetics  model.  The  approach  followed  here  consists  of  building  a  detailed  chemical 
kinetics  model  (Section  3.3),  followed  by  its  validation  with  relevant  experimental  data 
published  in  the  literature  (Section  3.4).  Once  validated,  the  model  is  then  used  to  interpret 
the  main  features  of  the  hypergolic  ignition  studied  in  this  work  through  various  sensitivity 
analyses  (Section  3.5). 

3.2.  Mixtures  under  consideration 

By  considering  the  equations  given  in  Ref.  1 4,  for  one  mole  of  DNB  decomposed  in  the  liquid 
phase,  there  will  be  an  equal  amount  of  HNCO,  C02,  H20  and  twice  this  amount  of  N20  in 
the  gas  phase.  N20  is  known  not  only  as  an  oxidant  but  also  as  a  monergol,  that  is,  its 
decomposition  is  so  exothermic  that  it  can  sustain  a  decomposition  flame.  The  amount  of 
gas-phase  HNO3,  at  a  given  temperature,  will  depend  on  the  type  of  liquid  solution  formed 
with  IL  and  WFNA;  either  an  ideal  solution,  in  which  case,  Raoult’s  law  will  be  followed,  or 
a  non-ideal  solution,  for  which,  the  values  of  the  activity  coefficients  of  the  species  would  be 
needed.  Furthermore,  here  the  reacting  liquid  mixture  is  inhomogeneous  in  nature. 
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Moreover,  a  small  amount  of  NO2  can  be  present  in  the  WFNA.  Therefore,  in  the  mixtures 
under  consideration  here,  there  is  one  type  of  fuel  molecule  (HNCO),  three  types  of  oxidizer 
molecules  (HNO3,  NO2  and  N20),  and  two  types  of  almost  unreactive  species  (H2O  and  CO2) 
prior  to  ignition.  Since  the  quantitative  composition  is  not  known,  various  mixtures  need  to 
be  considered  parametrically  in  the  kinetics  model  simulations  such  as:  1  HNCO/  1  CO2  /  1 
H20  /  2  N20  /  a  HN03  /  0.01  N02,  with  0.1<  a  <  2. 

3.3.  Building  of  a  detailed  chemical  kinetics  model 

Dagaut  et  al.  [281  have  published  a  detailed  kinetics  model  of  HCN  oxidation  (with  41  species 
and  250  reactions)  in  which  an  HNCO  sub-mechanism  has  been  considered.  Their  HNCO 
and  HONO2  sub-mechanisms  were  rather  limited,  therefore,  here  additional  reactions  have 
been  added  to  the  Dagaut  et  al.  model.  Also,  for  completeness,  the  C-species  sub-model  has 
been  included,  although  its  impact  to  the  ignition  modeling  of  HNCO  is  expected  to  be 
minimal  as  it  is  for  HCN  oxidation.  Thermal  decompositions  of  all  the  reactants  (HNCO, 
HNO3,  etc.)  are  also  explicitly  considered.  The  detailed  chemical  kinetics  model  written  for 
this  work  contains  70  species  and  506  reactions.  This  model  is  reported  as  supplementary 
information.  Rate  coefficients  employed  are  from  various  literature  sources. [28-30]  It  is  to 
be  noted  that  the  reaction  kinetics  of  some  of  the  steps  of  isocyanic  acid  (HNCO) 
decomposition  and  oxidation  in  the  gas  phase  have  previously  been  extensively  studied 
because  HNCO  is  the  decomposition  product  of  cyanuric  acid,  which  is  the  prime  reagent  of 
the  RapreNOx  process  for  NOx  removal.  [29]  However,  despite  these  studies,  no  detailed 
kinetics  model  has  been  specifically  devoted  to  HNCO  and  validated  with  all  of  the  available 
experimental  data  in  the  literature,  and  this  work  attempts  to  bridge  this  gap  in  our  model. 
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3.4.  Validation  of  the  detailed  chemical  kinetics  model 


The  CHEMKIN-II  [31]  and  SENKIN  [32]  codes  were  used  to  integrate  the  time-dependent 
rate  equations  derived  from  the  reaction  mechanism.  Direct  experimental  data  reported  in  the 
literature  concerning  various  HNCO-based  systems,  such  as  HNCO/Ar,  HNC0/N20/Ar,  and 
EINCO/NCVAr  are  considered  here.  These  data  were  not  considered  in  the  HNCO 
submechanism  by  Dagaut  et  al.  [28]  for  the  validation  of  their  HCN  oxidation  kinetics  model, 
which  was  instead  validated  with  available  data  on  HCN-based  systems  such  as 
HCN/EINCVAr  and  HCN/CVAr.  For  model  validation  purposes,  comparisons  between 
experimental  data  and  predictions  of  the  models  are  reported  in  Figs.  2-6  for  HCN-based 
systems.  As  shown  in  these  figures,  our  expanded  model  results  in  little  changes  to  Dagaut’ s 
predictions  for  HCN-based  systems. 
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Fig.  2:  Calculated  and  experimental  OH  profiles  for  the  mixture;  10.3  ppm  HNO3  + 
0.34%  HCN  in  Ar,  T  =  1492  K  and  P  =  1.01  atm.  Points  are  experimental  values  (as 
reported  in  Ref.  28),  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 
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Time  (|js) 

Fig.  3:  Calculated  and  experimental  N  profiles  for  the  mixture;  50  ppm  HCN  +  50  ppm 
O2  in  Ar,  T  =  2600  K  and  P  =  1.78  atm.  Points  are  experimental  values  (as  reported  in 
Ref.  28),  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 
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Fig.  4:  Calculated  and  experimental  O  profiles  for  the  mixture;  50  ppm  HCN  +  50  ppm 
O2  in  Ar,  T  =  2826  K  and  P  =  1.69  atm.  Points  are  experimental  values  (as  reported  in 
Ref.  28),  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 


Time  (|js) 

Fig.  5:  Calculated  and  experimental  N  profiles  for  the  mixture;  100  ppm  HCN  +  1000 
ppm  O2  in  Ar,  T  =  2718  K  and  P=1.78  atm.  Points  are  experimental  values  (as  reported 
in  Ref.  28),  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 
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Fig.  6:  Calculated  and  experimental  H  profiles  for  the  mixture;  50  ppm  HCN  +  250 
ppm  O2,  T  =  2625  K  and  P  =  1.82  atm.  Points  are  experimental  values  (as  reported  in 
Ref.  28),  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 

More  importantly,  this  expanded  model  considers  the  low  temperature  (673-773  K)  data  that 
are  available  for  HNCO/NCVAr  mixtures.  [30]  The  temperature  reached  in  the  gas  phase 
when  hypergolic  ignition  occurs  is  estimated  to  be  about  423  K  and  therefore  these  relatively 
low  temperature  data  are  particularly  important  for  the  validation  of  the  present  detailed 
chemical  kinetics  model  and  for  its  predictive  ability  at  around  423  K.  However,  this 
temperature  range  was  beyond  the  scope  of  Dagaut  et  al.  study  and  almost  all  of  the  data  they 
considered  were  obtained  at  900  K  or  above.  The  Dagaut  et  al.  model  does  not  correctly 
predict  the  HNCO/NO2  reactivity  in  the  673-773  K  temperature  range  (Figure  7),  whereas  this 
model  comes  much  closer  to  predicting  the  experimental  HNCO  mole  fraction  by  using  the 
rate  coefficients  recommended  by  He  et  al.  [30] 
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Fig.  7:  Calculated  and  experimental  HNCO  profiles  for  the  mixture;  1275.84  ppm 
HNCO  +  1485.68  ppm  NO2  in  Ar,  T  =  673  K  and  P  =  0.658  atm.  Points  are  experimental 
values  of  Ref.  30,  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 


The  predictions  of  the  present  model  are  compared  to  low  temperature  species  profiles 
obtained  in  a  static  reactor  (Figs.  8a-f).[30]  Experimental  data  and  predicted  data  agree  well 
for  all  the  species  profiles  except  for  CO2  whose  formation  is  generally  over  predicted. 
However,  experimental  data  and  predicted  data  remain  within  a  factor  2. 
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Fig.  8a:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 

NO  (A  and _ ),  C02  (❖  and _ ),  HNCO  (O  and _ ),  N20  (X  and 

and  N02  (fand . )  for  the  mixture;  1275.84  ppm  HNCO  +  1485.68 

ppm  N02  in  Ar,  T  =  673  K  and  P  =  0.658  atm. 


Fig.  8b:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 
C02  (O  and _ ),  HNCO  (O  and _ ),  N20  (X  and _ ),  and  N02 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


16 


(♦and . )  for  the  mixture;  1262.4  ppm  HNCO  +  1460.8  ppm  NO2  in  Ar,  T  =  723 

K  and  P  =  0.658  atm. 


Time  (s) 

Fig.  8c:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 

NO  (A  and _ ),  C02  (O  and _ ),  HNCO  (O  and _ ),  N20  (X  and 

and  N02  (iand . )  for  the  mixture;  1272.6  ppm  HNCO  +  1475  ppm 

N02  in  Ar,  T  =  773  K  and  P  =  0.658  atm. 
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Fig.  8d:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 


C02  (O  and _ ),  HNCO  (O  and _ ),  N20  (X  and _ ),  and  NOz 

(♦and . )  for  the  mixture;  2687  ppm  HNCO  +  1280  ppm  N02,  T  =  723  K  and  P  = 

0.658  atm. 


Fig.  8e:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 
C02  (O  and _ ),  HNCO  (O  and _ ),  N20  (X  and _ ),  and  N02 
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(♦and . )  for  the  mixture;  1316.5  ppm  HNCO  +  2858.4  ppm  NO2,  T  =  723  K  and 

P  =  0.658  atm. 


Fig.  8f:  Experimental  points  (from  Ref.  30)  and  calculated  curves  (present  model)  for 

NO  (A  and _ ),  C02  (O  and _ ),  HNCO  (O  and _ ),  N20  (X  and 

N02  (4and . )  for  the  mixture;  1343.5  ppm  HNCO  +  5274.8  ppm 

N02  in  Ar,  T  =  723  K  and  P  =  0.658  atm. 


The  present  study  is  devoted  to  hypergolic  ignition  chemistry  which  occurs  at  ambient  or 
relatively  low  temperatures.  However,  for  completeness,  high-temperature  data  above  1000 
K  have  also  been  considered.  The  agreement  between  experiments  and  calculations  is  good 
for  almost  all  of  the  experimental  data,  [28,  33-36]  whatever  the  system  studied,  and  whatever 
the  temperature  and  pressure  range  employed.  In  particular,  as  shown  in  Figs.  9a-9d,  the 
NCO  profiles  at  high  temperature  are  generally  better  predicted  with  the  present  model  than 
with  the  model  of  Dagaut  et  al.  [28]  This  is  very  important  because  the  amounts  of  NCO 
produced  are  generally  quantitatively  important  on  the  order  of  hundreds  of  ppm. 
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Fig.  9:  Calculated  and  exp enmental  NCO  profiles  for  the  mixtures:  a)  1  mol%  HNCO  +  1  mol%  N20  in  Ar,  T  =  2380  K  and  P  =  1.07 
atm.  b)  0.3  mol%  HNCO  + 1  mol%  N20  in  Ar,  T  =  2380  K  and  P  =  1. 14  atm.  c)  0.3  mol%  HNCO  +1.5  mol%  NO  in  Ar,  T  =  2383  K 
andP  =  1.12  atm.  cl)  0.2  mol%  HNCO  in  Ar,  T  =  2920  K  andP  =  0.92  atm.  Points  are  experimental  values  of  Ref.  33,  solid  line  is  from 
present  model  and  dashed-line  is  from  Ref.  28. 
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In  contrast,  the  NH  profiles  are  generally  better  predicted  when  using  the  model  of  Dagaut  et 
al.  (Fig.  10a)  but  it  is  not  always  the  case  (Fig.  10b).  Flowever,  the  amounts  of  NH  are  rather 
low,  10  ppm  at  most,  and  predictions  using  the  present  model  are  still  within  a  factor  of  2. 
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Fig.  10a:  Calculated  and  experimental  NH  profiles  for  the  mixture,  0.1  mol%  HNCO  + 
0.24  mol%  N2O  in  Ar,  T  =  2175  K  and  P  =  0.336  atm.  Points  are  experimental  values  of 
Ref.  35,  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 
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Fig.  10b:  Calculated  and  experimental  NH  profiles  for  the  mixture:  17  ppm  HNCO  in 
Ar,  T  =  3186  K  and  P  =  0.805  atm.  Points  are  experimental  values  of  Ref.  36,  solid  line  is 
from  present  model  and  dashed-line  is  from  Ref.  28. 


For  OH  radicals,  as  shown  in  Fig.  11,  both  models  predict  similar  profiles  and  are  in 
agreement  with  the  experimental  data. 
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Fig.  11:  Calculated  and  experimental  OH  profiles  for  the  mixture:  0.5  mol%  HNCO  + 
0.25  mol%  N2O  in  Ar,  T  =  2555  K  and  P  =  0.505  atm.  Points  are  experimental  values  of 
Ref.  35,  solid  line  is  from  present  model  and  dashed-line  is  from  Ref.  28. 


The  occasional  disagreement  observed  at  high  temperatures  suggests  further  refinements  and 
studies  are  needed,  and  reassessment  of  some  of  the  rate  coefficients  used. 


3. 5.  Predictions  of  the  kinetics  model 


The  full  kinetics  model  with  reaction  R506  shown  below  but  not  originally  considered,  was 
used  to  predict  the  fate  of  premixed  gas  mixtures;  1  HNCO/  1  CO2  /  1  H20  /  2  N20/  a  HNO3/ 
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0.01  N02  with  0.1<  a  <  2  at  initially  423  K  and  1  atm,  as  well  as  examine  the  ignition  delays. 
The  model  is  able  to  simulate  ignition  at  the  relatively  low  temperature  of  423  K.  However, 
computed  ignition  delay  times  are  longer  by  several  orders  of  magnitude  compared  to 
reported  experimental  values.  [14,16]  The  cause  of  this  is  unlikely  to  be  erroneous 
extrapolations  of  the  rate  coefficients  to  423  K  because  the  model  is  reasonably  predictive  in 
the  wide  673-3186  K  temperature  range.  Therefore,  additional  reactions  need  to  be  included 
in  the  model.  Wooldridge  et  al.  reported  that  isocyanic  acid  and  nitric  acid  react  at  room 
temperature  to  form  C02,  H20  and  N20.[37] 

H-N=C=0  +  H0-N02  -*■  C02  +  H20  +  N20  (R506) 

The  reaction  probably  proceeds  via  02N-N=C=0  formation  and  H20  elimination,  with  the 
former  rapidly  decomposing  to  give  C02  and  N20.  The  above  global  reaction  was  thus 
incorporated  into  the  kinetics  model  with  a  rate  coefficient  value  sufficiently  high  to  compute 
ignition  delay  times  that  were  comparable  to  experimentally  observed  values,  but  sufficiently 
low  to  have  no  influence  on  the  reported  predictions  of  Figs  2-11.  For  example,  ignition 
delay  times  x  of  ~700  ms  at  constant  volume  and  of  ~2  s  at  constant  pressure  for  a  =  1  were 
computed  (see  Table  1  and  Fig.  12).  No  additional  reactions  seem  to  explain  this  feature. 


a 

x  (s)  at  constant  volume 

x  (s)  at  constant  pressure 

0.05 

No  ignition 

No  ignition 

0.1 

No  ignition 

No  ignition 

0.2 

2500 

No  ignition 

0.3 

9 

275 

0.4 

2 

12 
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0.5 

1.2 
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0.6 

1 

3 

0.7 

0.8 

2 

0.8 

0.8 

2 

0.9 

0.75 

2 

1 

0.7 

2 

1.5 

0.8 

3 

Table  1.  Calculated  ignition  delays  x  at  constant  pressure,  and  at  constant  volume  for 
the  mixtures;  1  HNCO/ 1  C02  / 1  H20  /  2  N20/  a  HNO3/  0.01  N02  with  0.05  <  a  <  1.5 
initially  at  423  K  and  1  atm. 


In  Table  1,  from  a  =  1  to  a  =  1.5,  x  begins  to  increase  with  increasing  HNO3  concentration 
because  here  the  excess  HNO3  does  not  play  the  role  of  an  oxidizer  but  of  a  diluent  which 
inhibits  the  reaction,  not  chemically  but  physically  by  absorbing  part  of  the  exothermicity. 


Fig.  12:  Computed  temperature  profile  for  the  mixture, 

HNC0/N02/N20/C02/H20/HNC>3:  1/0.01/2/1/1/1.  The  system  is  assumed  to  be  adiabatic 
and  at  constant  volume.  P  =  1  atm  and  T  =  423  K.  Ignition  delay  x  is  730  millisecond. 
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3.6.  Summary  of  gas-phase  HNCO/HONO  3/N 2O/H2O/CO 2/NO 2  chemistry 


Sensitivity  analyses,  performed  with  the  full  kinetics  model,  enable  the  identification  of 
redundant  species  and  reactions,  that  is,  species  and  reactions  that  do  not  affect  the  flow  of  the 
reaction  system.  Several  sensitivity  analyses  exist.  In  the  method  of  brute  force  sensitivity 
analysis,  the  kinetics  code  is  run  by  perturbing  one  parameter  and  comparing  the  results  to  the 
original  unperturbed  model.  For  instance,  a  rate  coefficient  can  be  set  equal  to  zero  (as  if  to 
remove  the  reaction  from  the  system)  and  see  if  the  reaction  is  important  in  the  process.  In 
another  approach,  the  value  of  the  rate  coefficient  of  a  reaction  may  be  multiplied  or  divided 
generally  between  two  and  five  to  ascertain  its  importance  in  the  model.  Modification  of  the 
rate  coefficient  of  an  important  reaction  will  induce  relatively  large  variations  in  the  predicted 
ignition  delay.  Such  sensitivity  analyses  were  performed  for  the  representative  mixture  in 
Fig.  12  during  pre-ignition  from  time  zero  to  ~700  ms.  Together  with  these  sensitivity 
analyses,  local  sensitivity  analyses  (as  implemented  in  CHEMKIN-II),  for  several  species  not 
reported  here  and  for  the  temperature  were  performed.  The  normalized  sensitivity  coefficient, 
co j,  for  temperature  is  defined  as: 


f  dT  1  A. 

'  Oh  Jr 

with  j  as  the  index  of  the  reaction  and  A  as  the  pre-exponential  factor  of  its  reaction  rate 

coefficient.  Essentially,  co y  is  a  measure  of  how  the  system  temperature  depends  upon  Aj  , 
where  a  positive  co,  indicates  an  increase  in  system  temperature  with  an  increase  in  the 
reaction  pre-exponential  factor,  and  a  negative  co,  indicates  a  decrease  in  system  temperature 
with  increasing  reaction  pre-exponential  factor.  From  these  sensitivity  studies,  the  gas-phase 
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chemistry  can  be  summarized  as  in  Fig.  13.  Reaction  506  (hereafter  noted  as  R506)  shows 
the  most  sensitive  in  the  model  even  at  the  shortest  times  and  lowest  temperatures.  The 
exothermicity  of  R506  (H-N=C=0  +  HO-NO2  — »  CO2  +  FFO  +  N2O)  induces  significant 
thermal  decomposition  of  nitric  acid  according  to  -R210: 


HONO2  (+M)  -*■  NO2+OH  (+M)  (-R210) 


The  OH  and  NO2  chain  carriers  may  then  react  according  to: 


HNCO  +  OH  — s 

►  NCO  +  H20 

(R40) 

HONO2  +  OH  - 

» H20  +  NO3 

(R475) 

NO  +  OH  (+M) 

-►  HONO  (+M) 

(R164) 

N02  +  OH  — >  NO  +  H02 

(R165) 

HNO  +  N02  - 

■>  HONO  +  NO 

(R153) 

NCO  +  N02  - 

->  CO  +  2NO 

(R70) 

NCO  radicals,  mostly  formed  by  R40,  also  react  with  NO2  according  to: 


NCO  +  N02  -*■  C02  +  N20 


(R71) 
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This  reaction  is  highly  exothermic  (~500  kJ  mol"1)  and  could  be  expected  to  play  an  important 
role  in  the  ignition  process.  However,  this  expectation  is  not  realized  because  model  does  not 
appear  to  be  very  sensitive  to  it  shown  Fig.  13. 

NO,  formed  by  R153,  R70  and  R165,  reacts  with  NO3  according  to: 


NO  +  N03  -» N02  +  N02  (-R178) 


where  the  NO3  is  formed  in  R475  as  well  as  in  R505. 


HNNO  +  N02  -» NNH  +  N03  (R505) 

The  HNO  radicals  reacting  in  R153  are  mostly  formed  by: 

HNCO  +  N03  -*■  C02  +  NO  +  HNO  (R479) 

HNNO  +  N02  -*■  HNO  +  NO  +  NO  (-R490) 

HNNO  species  reacting  in  R505  and  -R490  are  formed  by: 

HNCO  +  N02  -*■  C02  +  HNNO  (R480) 

HO2  is  formed  by  reaction  R165  and  reacts  according  to  R175: 

N02  +  H02  ->  HONO  +  02  (R175) 
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H2O  and  N2O,  which  are  initially  present,  do  not  react  during  the  ignition  delay  period 
because  of  the  low  temperatures  but  are  respectively  formed  by  R506  and  R475,  and  by  R506 
and  R71. 
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Fig.  13:  Sensitivity  plot  for  temperature  for  the  five  most  sensitive  reactions:  (R506)  H- 
N=C=0  +  HO-NO2  CO2  +  H20  +  N20,  (R40)  HNCO  +  OH  — >  NCO  +  H20,  (R480) 

HNCO  +  N02  C02  +  HNNO,  (R165)  NO  +  H02  N02  +  OH  and  (R475)  HONO2  + 

OH  ->  H20  +  N03. 


Fig.  14  shows  the  model’s  fidelity  in  predicting  the  ignition  delay  time  when  R506  rate 
coefficient  value  is  doubled  and  when  it  is  halved.  The  predicted  ignition  delay  values  are  ~ 
350  ms  and  1400  ms.  Other  reactions  are  less  influential  on  the  predicted  ignition  delays. 
This  further  indicates  the  importance  of  R506. 
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Time  (s) 

Fig.  14:  Brute-force  sensitivity  analysis.  Comparison  between  calculated  ignition  delays 
with  variation  of  the  rate  coefficient  for  R506.  Solid  line:  present  model.  Dotted  line: 
pre-exponential  factor  multiplied  by  2.  Dashed  line:  pre-exponential  factor  divided  by  2. 

4.  Conclusions 

Experimentally  demonstrated  hypergolicity  between  an  ionic  liquid  (IL)  and  nitric  acid  has 
been  interpreted  here  by  using  a  detailed  chemical  kinetics  model.  An  exothermic  reaction 
between  isocyanic  acid  and  nitric  acid  needs  to  be  considered  in  the  model  to  explain  ignition. 
The  full  kinetics  model  is  shown  to  be  valid  in  a  wide  temperature  range  (673-3186  K). 
Various  sensitivity  analyses  show  that  this  exothermicity  induces  the  thermal  decomposition 
of  some  of  the  unreacted  nitric  acid  to  generate  OE1  and  NO2,  which  then  act  as  chain  carriers 
up  to  autoignition.  Several  other  reactions,  important  below  1000  K,  have  also  been  added  in 
the  present  model.  The  present  full  model  with  70  species  and  506  reactions  can  be  cast  into  a 
reduced  form  and  still  capture  accurately  the  overall  chemical  description  of  the  ignition 
process.  Reduced  forms  often  can  eliminate  ~  90%  of  the  original  reactions  from  the  full 
model  and  this  can  allow  for  the  implementation  of  chemistry  into  computational  codes  which 
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deal  with  the  physics  of  hypergolicity  (i.e.,  fluid  mechanical  and  thermal  transfer  effects).  It 


appears  from  the  results  of  the  present  study  that  the  published  chemical  kinetics  model  of 
HNCO  oxidation  needs  to  be  further  investigated,  and  that  reappraisal  of  some  of  the  rate 
coefficients  will  probably  be  necessary. 
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7.  Supplemental  Information 

Table  of  detailed  chemical  kinetics  reaction  model  used  in  this  work  to  simulate  ignition  delay 
times. 


(k  =  A.  Tb.  exp(-E/RT)) 


REACTIONS  CONSIDERED 

A 

b 

E 

1.  HCN+M  =  H+CN+M 

3.40E+35 

-5.1 

133000.0 

N2  Enhanced  by 

O.OOOE+OO 

02  Enhanced  by 

1.500E+00 

H20  Enhanced  by 

1.000E+01 

2.  HCN+N2=H+CN+N2 

3.60E+26 

-2.6 

124890.0 

3.  HCN+M=HNC+M 

1.60E+26 

-3.2 

54600.0 

AR  Enhanced  by 

7.000E-01 

H20  Enhanced  by 

7.000E+00 

C02  Enhanced  by 

2.000E+00 

4.  CN+H2=HCN+H 

1.10E+05  : 

2.6 

1908.0 

5.  HCN+0=NCO+H 

1.40E+04 

2.6 

4980.0 

6.  HCN+OCN+OH 

4.20E+10 

0.4 

20665.0 

7.  HCN+0=NH+CO 

3.50E+03 

2.6 

4980.0 

8.  HCN+OH=CN+H20 

3.90E+06 

1.8 

10300.0 

9.  HCN+OH=HOCN+H 

5.90E+04 

2.4 

12500.0 

10.  HCN+OH=HNCO+H 

2.00E-03 

4.0 

1000.0 

11.  HCN+OH=NH2+CO 

7.80E-04 

4.0 

4000.0 

12.  HCN+02=CN+H02 

3.00E+13 

0.0 

75100.0 
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13.  HCN+CN=NCCN+H 

14.  HN C+H=HCN +H 

15.  HNC+0=NH+C0 

16.  HNC+OH=HNCO+H 

17.  HNC+CN=NCCN+H 

18.  CN+0=C0+N 

19.  CN+OH=NCO+H 

20.  CN+02=NC0+0 

21.  CN+02=NC0+0 

22.  CN+02=N0+C0 

23.  CN+C02=NC0+C0 

24.  CN+CH20=HCN+HC0 

25.  CN+N0=NC0+N 

26.  CN+N02=NC0+N0 

27.  CN+N02=C0+N20 

28.  CN+N02=N2+C02 

29.  CN+HN0=HCN+N0 

30.  CN+H0N0=HCN+N02 

31.  CN+N20=NCN+N0 

32.  CN+HNC0=HCN+NC0 

33.  CN+NC0=NCN+C0 

34.  HNC0+M=C0+NH+M 

N2  Enhanced  by 

35.  HNCO+H=NH2+CO 

36.  HNCO+H=NCO+H2 


1.50E+07 

1.7 

1530.0 

7.80E+13 

0.0 

3600.0 

4.60E+12 

0.0 

2200.0 

2.80E+13 

0.0 

3700.0 

1.00E+13 

0.0 

0.0 

1.90E+12  0.5  723.0 

1.00E+15  -0.4  0.0 

7.20E+12  0.0  -417.0 

********  -2  0  00 

2.80E+17  -2.0  0.0 

3.70E+06  2.2  26900.0 
4.20E+13  0.0  0.0 

9.60E+13  0.0  42100.0 


5.30E+15  -0.8 

344.0 

4.90E+14  -0.8 

344.0 

3.70E+14  -0.8 

344.0 

1.80E+13  0.0 

0.0 

1.20E+13  0.0 

0.0 

3.80E+03  2.6 

3700.0 

1.00E+13  0.0 

0.0 

1.80E+13  0.0 

0.0 

1.10E+16  0.0  86000.0 


3.60E+04  2.5  2345.0 

9.00E+07  1.7  13900.0 


1.500E+00 
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37.  HNC0+0=NC0+0H 

2.20E+06 

2.1 

11430.0 

38.  HNC0+0=NH+C02 

9.60E+07 

1.4 

8520.0 

39.  HNC0+0=HN0+C0 

1.50E+08 

1.6 

44012.0 

40.  HNC0+0H=NC0+H20 

3.60E+07 

1.5 

3600.0 

41.  HNC0+H02=NC0+H202 

3.00E+11  0.0  22000.1 

42.  HNC0+02=HN0+C02 

1.00E+12 

0.0 

35000.0 

43.  HNCO+NH=NH2+NCO 

3.00E+13 

0.0 

23700.0 

44.  HOCN+H=HNCO+H 

3.10E+08 

0.8 

1917.0 

45.  HOCN+H=NH2+CO 

1.20E+08 

0.6 

2076.0 

46.  HOCN+H=H2+NCO 

2.40E+08 

1.5 

6617.0 

47.  HOCN+OOH+NCO 

1.70E+08 

1.5 

4133.0 

48.  H0CN+0H=H20+NC0 

1.20E+06 

2.0 

-248.0 

49.  HOCN+NH2=NCO+NH3 

9.20E+05 

i  1.9 

3646.0 

50.  HCN0=HCN+0 

4.20E+31  -6.1  61210.0 

51.  HCNO+H=HCN+OH 

7.20E+10 

0.8 

8612.0 

52.  HCN0+0=HC0+N0 

6.30E+13 

0.0 

0.0 

53.  HCN0+0H=CH20+N0 

1.00E+12 

0.0 

0.0 

54.  HCN0+0=NC0+0H 

7.00E+12 

0.0 

0.0 

55.  HCN O+O H=N O+C 0+H2 

6.50E+12 

0.0 

0.0 

56.  HCNO+OH=NCO+H+OH 

4.50E+12 

0.0 

0.0 

57.  HCNO+OH=NCO+H20 

3.50E+12 

0.0 

0.0 

58.  HCNO+OH=HCO+HNO 

4.50E+12 

0.0 

0.0 

59.  NCO+M=N+CO+M 

2.20E+14  0.0  54050.0 

N2  Enhanced  by  1.500E+00 

60.  NCO+H=CO+NH 

7.20E+13  0.0 

1000.0 
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0.0 


61.  NC0+0=N0+C0 

62.  NCO+OH=HON+CO 

63.  NCO+OH=H+CO+NO 

64.  NC0+H02=HNC0+02 

65.  NC0+02=N0+C02 

66.  NCO+HCO=HNCO+CO 

67.  CH20+NC0=HNC0+HC0 

68.  NC0+N0=N20+C0 

69.  NC0+N0=N2+C02 

70.  NC0+N02=C0+N0+N0 

71.  NC0+N02=C02+N20 

72.  NCO+HNO=HNCO+NO 

73.  NC0+H0N0=HNC0+N02 

74.  NCO+N=N2+CO 

75.  NCO+NH3=HNCO+NH2 

76.  NCO+NCO=CO+CO+N2 

77.  H2CN=HCN+H 

78.  H2CN+H=HCN+H2 

79.  H2CN+0=HCN+0H 

80.  H2CN+0H=HCN+H20 

81.  H2CN+0H=HCN+H20 

82.  H2CN+02=CH20+N0 

83.  H2CN+NH2=HCN+NH3 

84.  HCNH=HCN+H 

85.  HCNH+H=H2 CN+H 


2.00E+15  -0.5 
5.30E+12  -0.1  5126.0 

8.30E+12  -0.1  18042.0 

2.00E+13  0.0  0.0 

2.00E+12  0.0  20000.0 
3.60E+13  0.0  0.0 

6.00E+12  0.0  0.0 

4.00E+19  -2.2  1743.0 

1.50E+21  -2.7  1824.0 

2.50E+1 1  0.0  -707.0 

3.00E+12  0.0  -707.0 

1.80E+13  0.0  0.0 

3.60E+12  0.0  0.0 

2.00E+13  0.0  0.0 

2.80E+04  2.5  980.0 

1.80E+13  0.0  0.0 

4.00E+28  -6.0  29897.0 
2.40E+08  1.5  -894.0 

1.70E+08  1.5  -894.0 

1.50E+19  -2.2  2166.0 

1.20E+06  2.0  -1192.0 
3.00E+12  0.0  5961.0 

9.20E+05  1.9  -1152.0 

6.10E+28  -5.7  24271.0 
2.00E+13  0.0  0.0 
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86.  HCNH+H=HCN+H2 

87.  HCNH+0=HNC0+H 

88.  HCNH+0=HCN+0H 

89.  HCNH+0H=HCN+H2 0 

90.  NH3+M=NH2+H+M 

91.  NH3+H=NH2+H2 

92.  NH3+0=NH2+0H 

93.  NH3+0H=NH2+H20 

94.  NH3+H02=NH2+H202 

95.  NH2+H=NH+H2 

96.  NH2+0=HN O+H 

97.  NH2+0=NH+0H 

98.  NH2+0H=NH+H20 

99.  NH2+H02=H2N0+0H 

100.  NH2+H02=NH3+02 

101.  NH2+02=H2N0+0 

102.  NH2+02=HN0+0H 

103.  NH2+NH2=NH3+NH 

104.  NH2+NH=NH3+N 

105.  NH2+N=N2+H+H 

106.  NH2+N0=N2+H20 

107.  NH2+NO=NNH+OH 

108.  NH2+N02=N20+H20 

109.  NH2+N02=H2N0+N0 

110.  NH2+HNO=NH3+NO 


2.40E+08 

1.5 

-894.0 

7.00E+13 

0.0 

0.0 

1.70E+08 

1.5 

-894.0 

1.20E+06  2.0  -1192.0 

2.20E+16 

0.0 

93470.0 

6.40E+05  2.4  10171.0 

9.40E+06 

1.9 

6460.0 

2.00E+06 

2.0 

566.0 

3.00E+11 

0.0 

22000.0 

4.00E+13  0.0  3650.0 

6.60E+14  -0.5 

0.0 

6.80E+12  0.0 

0.0 

4.00E+06 

2.0 

1000.0 

5.00E+13 

0.0 

0.0 

9.20E+05 

1.9 

-1152.0 

2.50E+11 

0.5 

29586.0 

6.20E+07 

1.2 

35100.0 

5.00E+13 

0.0 

10000.0 

9.20E+05 

1.9 

2444.0 

7.00E+13 

0.0 

0.0 

2.80E+20 

-2.7 

1258.0 

2.30E+10 

0.4 

-814.0 

1.60E+16  -1.4  268.0 

6.50E+16  -1.4  268.0 

3.60E+06  1.6  -1250.0 
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111.  NH2+H0N0=NH3+N02 

112.  NH+H=N+H2 

113.  NH+0=N0+H 

114.  NH+OH=HN 0+H 

115.  NH+0H=N+H20 

116.  NH+02=HN 0+0 

117.  NH+02=N 0+0H 

118.  NH+NH=N2+H+H 

119.  NH+N=N2+H 

120.  NH+N0=N20+H 

121.  NH+N0=N20+H 

122.  NH+NON2+OH 

123.  NH+N02=N20+0H 

124.  NH+HON 0=NH2+N 02 

125.  N+0H=N0+H 

126.  N+02=N0+0 

127.  N+N0=N2+0 

128.  H2N0+M=HN0+H+M 

H20  Enhanced  by 

129.  H2NO+M=HNOH+M 

H20  Enhanced  by 

130.  H2NO+H=HNO+H2 

131.  H2NO+H=NH2+OH 

132.  H2N0+0=HN0+0H 

133.  H2N0+0H=HN0+H20 


7.10E+01  3.0  -4940.0 


3.00E+13  0.0 

0.0 

9.20E+13 

0.0 

0.0 

2.00E+13 

0.0 

0.0 

5.00E+11 

0.5 

2000.0 

4.60E+05 

2.0 

6500.0 

1.30E+06 

1.5 

100.0 

2.50E+13 

0.0 

0.0 

3.00E+13  0.0  0.0 

2.90E+14  -0.4  0.0 

-0.2  0.0 

-0.2  0.0 

0.0 
0.0 

0.0 

6280.0 
0.0 

-2.8  64915.0 


1.10E+29 

-4.0 

44000.0 

3.00E+07 

2.0 

2000.0 

5.00E+13 

0.0 

0.0 

3.00E+07 

2.0 

2000.0 

2.00E+07  2.0  1000.0 


2.20E+13 

1.00E+13  0.0 
1.00E+13  0.0 
3.80E+13  0.0 
6.40E+09  1.0 

3.30E+12  0.3 
2.80E+24 

1.000E+01 

1.000E+01 
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134.  H2N0+H02=HN0+H202 


2.90E+04  2.7  1600.0 


135.  H2N0+02=HN0+H02 

3.00E+12  0.0 

25000.0 

136.  H2NO+NO=HNO+HNO 

2.00E+04  2.0 

13000.0 

137.  H2NO+NH2=HNO+NH3 

3.00E+12  0.0 

1000.0 

138.  H2NO+N 02=H0N O+HN 0 

6.00E+11  0.0  2000.0 

139.  HN OH+M=HN O+H+M 

2.00E+24  -2.8 

58934.0 

H20  Enhanced  by  1.000E+01 

140.  HNOH+H=NH2+OH 

4.00E+13  0.0 

0.0 

141.  HN OH+H=HN 0+H2 

4.80E+08  1.5 

378.0 

142.  HN 0H+0=HN O+OH 

7.00E+13  0.0 

0.0 

143.  HNOH+0=HNO+OH 

3.30E+08  1.5 

-358.0 

144.  HNOH+OH=HNO+H20 

2.40E+06  2.0 

-1192.0 

145.  HN OH+H02=HN 0+H2  02 

2.90E+04  2.7  -1600.0 

146.  HN OH+02=HN 0+H02 

3.00E+12  0.0 

25000.0 

147.  HNOH+NH2=NH3+HNO 

1.80E+06  1.9 

-1152.0 

148.  HN OH+N 02=HON O+HN O 

6.00E+11  0.0  2000.0 

149.  HNO+H=NO+H2 

4.40E+11  0.7  650.0 

150.  HNO+ONO+OH 

2.30E+13  0.0 

0.0 

151.  HNO+OH=NO+H20 

3.60E+13  0.0 

0.0 

152.  HN0+02=H02+N0 

2.00E+13  0.0  16000.0 

153.  HNO+N02=HONO+NO 

6.00E+1 1  0.0 

2000.0 

154.  HN0+HN0=N20+H20 

9.00E+08  0.0 

3100.0 

155.  HCO+HNO=NO+CH20 

6.00E+1 1  0.0 

2000.0 

156.  HON +M=N O+H+M 

5.10E+19  -1.7 

16045.0 

AR  Enhanced  by  7.000E-01 
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H20  Enhanced  by  7.000E+00 

C02  Enhanced  by  2.000E+00 


157.  HON +H=HN O+H 

2.00E+13 

0.0 

0.0 

158.  HON+H=OH+NH 

2.00E+13 

0.0 

0.0 

159.  HON +0=OH+N O 

7.00E+13 

0.0 

0.0 

160.  HON+OH=HONO+H 

4.00E+13 

0.0 

0.0 

161.  H0N+02=N02+0H 

1.00E+12 

0.0 

4968.0 

162.  H+NO(+M)=HNO(+M)  1.50E+15  -0.4  0.0 

Low  pressure  limit:  0.23000E+15  0.20600E+00  -0.15500E+04 


N2 

Enhanced  by 

1.000E+00 

H20 

Enhanced  by 

1.000E+01 

02 

Enhanced  by 

1.500E+00 

H2 

Enhanced  by 

2.000E+00 

C02 

Enhanced  by 

3.000E+00 

163.  N0+0(+M)=N02(+M)  1.30E+15  -0.8  0.0 

Low  pressure  limit:  0.75000E+20 -0.14100E+01  O.OOOOOE+OO 


N2 

Enhanced  by 

1.700E+00 

02 

Enhanced  by 

1.500E+00 

H20 

Enhanced  by 

1.000E+01 

164.  NO+OH(+M)=HONO(+M)  2.00E+12  -0.1  -721.0 

Low  pressure  limit:  0.50000E+24  -0.25100E+01  -0.68000E+02 


N2 

Enhanced  by 

1.000E+00 

02 

Enhanced  by 

1.000E+00 

H20 

Enhanced  by 

6.700E+00 

AR 

Enhanced  by 

6.700E-01 
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165.  N0+H02=N02+0H 

2.10E+12 

0.0 

-480.0 

166.  HCO+NOHNO+CO 

7.00E+13 

-0.4 

0.0 

167.  HON 0+H=HN O+OH 

5.60E+10 

0.9 

5000.0 

168.  HONO+H=NO+H20 

8.10E+06 

1.9 

3850.0 

169.  H0N0+0=N02+0H 

1.20E+13 

0.0 

6000.0 

4.00E+12  0.0  0.0 

3.50E-01  3.6  12100.0 

1.30E+14  0.0  362.0 

3.90E+12  0.0  -238.0 

1.30E+13  0.0  0.0 


170.  H0N0+0H=N02+H20 

171.  H0N0+H0N0=N0+N02+H20 

172.  N02+H=N0+0H 

173.  N02+0=N0+02 

174.  N02+0(+M)=N03(+M) 

Low  pressure  limit:  0.10000E+29 -0.40800E+01  0.24700E+04 
N2  Enhanced  by  1.500E+00 

02  Enhanced  by  1.500E+00 

H20  Enhanced  by  1.000E+01 

175.  N02+H02=H0N0+02 


176.  N02+H2=H0N0+H 

177.  N02+N02=N0+N0+02 

178.  N02+N02=N03+N0 

179.  C0+N02=C02+N0 

180.  HC0+N02=C0+N0+0H 

181.  HC0+N02=H+C02+N0 

182.  N03+H=N02+0H 

183.  N03+0=N02+02 

184.  N03+0H=N02+H02 

185.  N03+H02=N02+02+0H 


6.30E+08  1.2  5000.0 

4.50E+12  0.0  27600.0 
1.60E+12  0.0  26123.0 
9.60E+09  0.7  20900.0 
9.00E+13  0.0  33800.0 
1.20E+23  -3.3  2355.0 

8.40E+15  -0.8  1930.0 

6.00E+13  0.0  0.0 

l.OOE+13  0.0  0.0 

1.40E+13  0.0  0.0 

1.50E+12  0.0  0.0 
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186.  N03+N02=N0+N02+02 


5.00E+10  0.0  2940.0 


187.  NNH=N2+H  6.50E+07  0.0  0.0 


188.  NNH+H=N2+H2 

1.00E+14  0.0 

0.0 

189.  NNH+0=N20+H 

1.00E+14 

0.0 

0.0 

190.  NNH+0=N2+OH 

8.00E+13 

0.0 

0.0 

191.  NNH+0=NH+N  O 

5.00E+13 

0.0 

0.0 

192.  NNH+0H=N2+H20 

5.00E+13 

0.0 

0.0 

193.  NNH+02=N2+H02 

2.00E+14 

0.0 

0.0 

194.  NNH+02=N2+H+02 

5.00E+13 

0.0 

0.0 

195.  NNH+NH=N 2+NH2 

5.00E+13 

0.0 

0.0 

196.  NNH+NH2=N2+NH3 

5.00E+13 

0.0 

0.0 

197.  NNH+NO=N2+HNO 

5.00E+13 

0.0 

0.0 

198.  N20(+M)=N2+0(+M) 

1.30E+12 

0.0 

62570.0 

Low  pressure  limit:  0.40000E+15  0.00000E+00  0.56600E+05 


N2 

Enhanced  by 

1.700E+00 

02 

Enhanced  by 

1.400E+00 

C02 

Enhanced  by 

3.000E+00 

H20 

Enhanced  by 

1.200E+01 

199.  N20+H=N2+0H 

200.  N20+H=N2+0H 

201.  N20+0=N0+N0 

202.  N20+0=N2+02 

203.  N20+0H=N2+H02 

204.  N20+0H=HN0+N0 

205.  N20+N0=N02+N2 


3.30E+10 

0.0 

4729.0 

4.40E+14 

0.0 

19254.0 

9.20E+13 

0.0 

27679.0 

3.70E+12 

0.0 

15936.0 

1.30E-02  4.7  36560.0 
1.20E-04  4.3  25080.0 
5.30E+05  2.2  46280.0 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


43 


206.  C0+N20=N2+C02 


2.70E+11  0.0  20237.0 


207.  NCCN+M=CN+CN+M  1.10E+34  -4.3  130079.0 


N2 

Enhanced  by 

1.500E+00 

02 

Enhanced  by 

1.500E+00 

H2 

Enhanced  by 

1.500E+00 

H20 

Enhanced  by 

1.000E+01 

C02 

Enhanced  by 

3.000E+00 

208.  NCCN+0=CN+NC0  4.60E+12  0.0  8880.0 

209.  NCCN+OH=CN+HOCN  1.90E+11  0.0  2900.0 

210.  N02+0H(+M)=H0N02(+M)  2.40E+13  0.0  0.0 

Low  pressure  limit:  0.64000E+33  -0.54900E+01  0.23510E+04 


N2 

Enhanced  by 

1.000E+00 

AR 

Enhanced  by 

7.000E-01 

H20 

Enhanced  by 

6.000E+00 

211.  H+H+M= 

=H2+M 

7.31E+17  -1.0 

0.0 

H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 

C02 

Enhanced  by 

3.750E+00 

212.  0+0+M= 

=02+M 

1.14E+17  -1.0 

0.0 

H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 

C02 

Enhanced  by 

3.750E+00 

213.  0+H+M= 

=OH+M 

6.20E+16  -0.6 

0.0 

H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 
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C02 


Enhanced  by  3.750E+00 


214.  H2+02=0H+0H  1.70E+13  0.0  47780.0 

215. 0+H2=0H+H  3.87E+04  2.7  6260.0 

216.  H+02=0H+0  4.40E+14  -0.1  16812.0 

217.  H+02+M=H02+M  8.00E+17  -0.8  0.0 

H20  Enhanced  by  1.000E+01 

218.  H+0H+M=H20+M  8.62E+21  -2.0  0.0 


H20  Enhanced  by  1.625E+01 

CO  Enhanced  by  1.875E+00 

C02  Enhanced  by  3.750E+00 


219.  H2+0H=H20+H 

2.16E+08 

1.5 

3430.0 

220.  H20+0=0H+0H 

1.50E+10 

1.1 

17260.0 

221.  H02+0H=H20+02 

2.89E+13 

0.0 

-497.0 

222.  H02+0=0H+02 

1.81E+13 

0.0 

-400.0 

223.  H+H02=H2+02 

4.28E+13 

0.0 

1411.0 

224.  H+H02=OH+OH 

1.69E+14 

0.0 

874.0 

225.  H+H02=H2 O+O 

3.01E+13 

0.0 

1721.0 

226.  H02+H02=H202+02  4.08E+02  3.3  1979.0 

227.  OH+OH(+M)=H2 02 (+M)  7.22E+13  -0.4  0.0 

Low  pressure  limit:  0.221 10E+20 -0.76000E+00  0.00000E+00 
TROE  centering:  0.50000E+00  0.10000E+09  0.10000E-05 


228.  H202+0H=H02+H20 

5.80E+14 

0.0 

9557.0 

229.  H202+H=H02+H2 

1.70E+12 

0.0 

3750.0 

230.  H2 02+H=H2 O+OH 

1.00E+13 

0.0 

3590.0 

231.  H202+0=H02+0H 

2.80E+13 

0.0 

6400.0 
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232.  C0+H02=C02+0H 


1.50E+14  0.0  23650.0 


233.  C0+0H=C02+H  4.40E+06  1.5  -740.0 

234.  C0+0+M=C02+M  2.83E+13  0.0  -4540.0 


H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 

C02 

Enhanced  by 

3.750E+00 

235.  C0+02=C02+0 

2.53E+12  0.0  47700.0 

236.  HC 0+M=H+C O+M 

1.85E+17  -1.0  17000.0 

H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 

C02 

Enhanced  by 

3.750E+00 

237.  HCO+OH=CO+H20 

1.00E+14  0.0 

0.0 

238.  HCO+0=CO+OH 

3.00E+13  0.0 

0.0 

239.  HC0+0=C02+H 

3.00E+13  0.0 

0.0 

240.  HCO+H=CO+H2 

7.22E+13  0.0 

0.0 

241.  HC0+02=C0+H02 

4.72E+12  0.0 

250.0 

242.  HC0+H02=C02+0H+H 

3.00E+13  0.0 

0.0 

243.  HCO+HCO=CH20+CO 

1.80E+13  0.0 

0.0 

244.  HCO+HCO=H2+CO+CO 

3.00E+12  0.0 

i  0.0 

245.  CH20+M=HCO+H+M 

1.26E+16  0.0 

77898.0 

H20 

Enhanced  by 

1.625E+01 

CO 

Enhanced  by 

1.875E+00 

C02 

Enhanced  by 

3.750E+00 

246.  CH20+H02=HC0+H202 

4.00E+12 

0.0 

11665.0 

247.  CH20+0H=HC0+H20 

1.72E+09 

1.2 

-447.0 
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248.  CH20+0=HC0+0H 


1.81E+13  0.0  3088.0 


249.  CH20+H=HC0+H2  1.26E+08  1.6  2170.0 

250.  CH20+02=HC0+H02  2.04E+13  0.0  39000.0 


251.  HNC0+0H=NH2+C02 

6.30E+10 

-0.1 

11643.0 

252.  NC0+0=N+C02 

8.00E+12  0.0 

2503.0 

253.  HNO+CH3=NO+CH4 

8.20E+05 

1.9 

953.0 

254.  NO+C=CO+N 

1.70E+13  0.0 

0.0 

255.  NO+C=CN+0 

1.10E+13  0.0 

0.0 

256.  NH2+N02=NH2N02 

3.50E+31  -6.8  3726.0 

257.  CH30+N0=CH20+HN0 

8.40E+12 

0.0 

2050.0 

258.  CH3+H02=CH30+0H 

2.28E+13 

0.0 

0.0 

259.  CH3+N02=CH30+N0 

1.30E+13 

0.0 

0.0 

260.  N2 04(+M)=N 02+N 02(+M) 

4.05E+18 

-1.1 

12840.0 

Low  pressure  limit:  0.19600E+29 -0.38000E+01  0.12840E+05 

261.  N02+H0N0=N0+H0N02  6.03E+01  0.0  0.0 

262.  N204+H20=H0N0+H0N02  2.52E+14  0.0  11586.0 

263.  C0+0H(+M)=C02+H(+M)  2.45E-03  3.7  -1242.8 

Low  pressure  limit:  0.11700E+08  0.13500E+01  -0.71700E+03 

264.  H0C0(+M)=C02+H(+M)  1.74E+12  0.3  32935.0 

Low  pressure  limit:  0.22900E+27 -0.30200E+01  0.35062E+05 

265.  HOCO(+M)=CO+OH(+M)  5.89E+12  0.5  33986.6 

Low  pressure  limit:  0.21900E+24 -0.18900E+01  0.35277E+05 

266.  CHCO+M=CO+CH+M  6.00E+15  0.0  58819.3 

267.  CH+0=C0+H  1.00E+14  0.0  0.0 

268.  CH+OH=HCO+H  3.00E+13  0.0  0.0 

Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


269.  CH+02=C0+0H 

270.  CH2S+0=C0+H+H 

271.  CH2S+0=C0+H2 

272.  CH2S+0H=CH20+H 

273.  CH2S+02=C0+0H+H 

274.  CH2T+0=C0+H+H 

275.  CH2T+0H=CH+H20 

276.  CH2T+0H=CH20+H 

277.  CH2T+02=C02+H+H 

278.  CH2T+02=CH2 0+0 

279.  CH2T+02=C02+H2 

280.  CH2T+02=C0+0H+H 

281.  CH2T+02=HC0+0H 

282.  CH3+0=CH20+H 

283.  CH3+0H=CH20+H2 

284.  CH3+0H=CH2T+H20 

285.  CH3+02=CH20+0H 

286.  CH3+02=CH302 

287.  CH3+HCO=CH4+CO 

288.  CH30H=CH3+0H 

289.  CH4+H02=CH3+H202 

290.  CH3+H02=CH4+02 

291.  CH4+0=CH3+0H 

292.  CH4+0H=CH3+H20 

293.  CH3CHOHCO+CH3 


2.00E+13  0.0  0.0 


1.50E+13 

0.0 

0.0 

1.50E+13  0.0 

0.0 

3.00E+13 

0.0 

0.0 

3.10E+13 

0.0 

0.0 

5.00E+13 

0.0 

0.0 

1.13E+07 

2.0 

3011.5 

2.50E+13 

0.0 

0.0 

1.60E+12 

0.0 

1003.8 

5.00E+13 

0.0 

9010.5 

6.90E+1 1 

0.0 

501.9 

8.60E+10 

0.0 

-501.9 

4.30E+10 

0.0 

-501.9 

8.43E+13  0.0 

0.0 

8.00E+12 

0.0 

0.0 

1.13E+06 

2.1 

2437.8 

5.20E+13 

0.0 

34894.8 

1.70E+60  -15.1 

18785.8 

3.20E+11 

0.5 

0.0 

2.80E+32  - 

5.0 

97992.3 

4.48E+13 

0.0 

24641.5 

1.00E+12 

0.0 

0.0 

1.90E+09  : 

1.4 

8675.9 

1.50E+06 

2.1 

2437.8 

2.00E+15 

0.0 

79110.9 
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294.  HC0+CH3 0H=CH2 0+CH20H 


9.60E+03  2.9  13097.5 


295.  CH20+CH3=HC0+CH4 

8.91E-13 

7.4 

-956.0 

296.  CH20+M=C0+H2+M 

8.20E+15 

0.0 

69550.7 

H2 

Enhanced  by 

2.900E+00 

02 

Enhanced  by 

1.200E+00 

N2 

Enhanced  by 

1.200E+00 

H20 

Enhanced  by 

1.850E+01 

CO 

Enhanced  by 

2.100E+00 

C02 

Enhanced  by 

4.300E+00 

297.  CH3 0+M=CH20+H+M  1.96E+37  -6.7  33269.6 


H2 

Enhanced  by 

2.900E+00 

02 

Enhanced  by 

1.200E+00 

N2 

Enhanced  by 

1.200E+00 

H20 

Enhanced  by 

1.850E+01 

CO 

Enhanced  by 

2.100E+00 

C02 

Enhanced  by 

4.300E+00 

298.  CH20H+M=CH20+H+M  1.54E+13  0.0  33054.5 


H2 

Enhanced  by 

2.900E+00 

02 

Enhanced  by 

1.200E+00 

N2 

Enhanced  by 

1.200E+00 

H20 

Enhanced  by 

1.850E+01 

CO 

Enhanced  by 

2.100E+00 

C02 

Enhanced  by 

4.300E+00 

299.  CH20H+HCH3+0H 

9.64E+13 

0.0 

0.0 

300.  CH20H+H=CH20+H2 

3.00E+13 

0.0 

0.0 

Distribution  A  :  Approved  for  public  release;  distribution  unlimited 
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301.  CH20H+0=CH20+0H 

302.  CH20H+0H=CH20+H20 

303.  CH20H+H02=CH20+H202 

304.  CH20H+CH3=CH20+CH4 

305.  CH20H+HC0=CH2  0+CH20 

306.  CH20H+HC0=CH30H+C0 

307.  CH20H+02=CH20+H02 

308.  CH20H+H202=CH3 0H+H02 

309.  CH30H+02=CH20H+H02 

310.  CH3+0H=CH30+H 

311.  CH3+02=CH30+0 

312.  CH3 0+H=CH20+H2 

313.  CH30+0=CH20+0H 

314.  CH30+0H=CH20+H20 

315.  CH30+02=CH20+H02 

316.  CH30+H02=CH20+H202 

317.  CH30+CH3=CH4+CH20 

318.  CH3 0+CH20=HC0+CH3 OH 

319.  CH302H=CH30+0H 

320.  CH3 OH+OH=C H2 OH+H2 O 

321.  CH30H+0H=CH30+H20 

322.  CH30H+H=CH20H+H2 

323.  CH30H+H=CH30+H2 

324.  CH3 0H+H=CH3+H20 

325.  CH3 0H+0=CH20H+0H 


5.00E+12  0.0  0.0 


5.00E+12 

0.0 

0.0 

1.20E+13 

0.0 

0.0 

2.40E+12 

0.0 

0.0 

1.80E+14 

0.0 

0.0 

1.20E+14 

0.0 

0.0 

1.80E+13  0.0  788.7 

3.01E+09  0.0  2581.3 

2.00E+13  0.0  44909.2 


5.74E+12 

-0.2 

13934.0 

4.30E+13 

0.0 

30807.8 

2.00E+13 

0.0 

0.0 

5.00E+12 

0.0 

0.0 

5.00E+12  0.0  0.0 

4.28E-13  7.6  -3537.3 
1.20E+13  0.0  0.0 

2.41E+13  0.0  0.0 

1.15E+11  0.0  1290.6 

6.46E+14  0.0  42997.1 
1.77E+04  2.6  -884.3 

1.77E+04  2.6  -884.3 


3.20E+13 

0.0 

6094.6 

8.00E+12 

0.0 

6094.6 

5.25E+12 

0.0 

5329.8 

3.88E+05 

2.5 

3083.2 

Distribution  A  :  Approved  for  public  release;  distribution  unlimited 
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326.  CH30H+CH30=CH20H+CH30H 

327.  CH30H+C2H3=CH20H+C2H4 

328.  CH302+H=CH30+0H 

329.  CH302+0=CH30+02 

330.  CH302+0H=CH30H+02 

331.  CH3 02+H02=CH3 02H+02 

332.  CH3 02+HC0=CH3 0+H+C02 

333.  CH302+CH3=CH30+CH30 

334.  CH3 02+H2 02=CH3 02H+H02 

335.  CH302+CH20=CH302H+HC0 

336.  CH302+CH302=CH30+CH30+02 

337.  CH302+CH302=CH30H+CH20+02 

338.  CH3 02+CH4=CH3 02H+CH3 

339.  CH3 02+H2=CH3 02H+H 

340.  CH302+CH30H=CH302H+CH20H 

341.  CH302H+0H=CH302+H20 

342.  CH302H+0=CH302+0H 

343.  C2H+0=C0+CH 

344.  C2H+OH=CHCO+H 

345.  C2H+02=CHC0+0 

346.  C2H+02=CH+C02 

347.  C2H+02=HC0+C0 

348.  C2H2+OH=C2H+H2  O 

349.  CHCO+H=CH2S+CO 

350.  CHC0+0=C0+C0+H 


1.51E+12  0.0  7002.9 

3.20E+01  3.2  7170.2 

9.60E+13  0.0  0.0 

3.60E+13  0.0  0.0 

6.00E+13  0.0  0.0 

4.60E+10  0.0  -2581.3 
3.00E+13  0.0  0.0 

2.00E+13  0.0  0.0 

2.40E+12  0.0  9942.6 

2.00E+12  0.0  11663.5 

7.80E+10  0.0  0.0 

1.30E+1 1  0.0  0.0 

1.80E+11  0.0  18475.1 

3.00E+13  0.0  26027.7 

1.80E+11  0.0  13719.0 

7.00E+12  0.0  1434.0 

2.80E+13  0.0  6405.3 

1.00E+13  0.0  0.0 


2.00E+13 

0.0 

0.0 

6.02E+11 

0.0 

0.0 

4.50E+15 

0.0 

25095.6 

2.41E+12 

0.0 

0.0 

2.71E+13 

0.0 

10492.3 

1.50E+14 

0.0 

0.0 

1.00E+14 

0.0 

0.0 

Distribution  A  :  Approved  for  public  release;  distribution  unlimited 
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351.  C2H2+0=CH2T+C0 

352.  C2H2+0=CHC0+H 

353.  C2H2+OH=HCCOH+H 

354.  C2H2+OH=CH2CO+H 

355.  C2H2+OH=CH3+CO 

356.  HCC0H+H=CH2C0+H 

357.  CH2CO+H=CHCO+H2 

358.  CH2C0+0=HC0+HC0 

359.  CH2C0+0H=CHC0+H20 

360.  CH2CO=CH2T+CO 

361.  C2H3+0H=C2H2+H20 

362.  C2H3+0=CH2C0+H 

363.  C2H3+02=CH20+HC0 

364.  CH3CO+H=CH2CO+H2 

365.  CH3C0+0=CH3+C02 

366.  CH3C0+CH302=CH3+C02+CH30 

367.  C2H4+H02=C2H3+H202 

368.  C2H4+02=C2H3  +H02 

369.  C2H3+CH30H=C2H4+CH30 

370.  C2H4+OCH3+HCO 

371.  C2H4+0=CH3C0+H 

372.  C2H4+0=C2H3+0H 

373.  C2H4+0H=C2H3+H20 

374.  CH3CHO+H=CH3CO+H2 

375.  CH3CH0+0=CH3C0+0H 


7.81E+03 

2.8 

501.9 

1.39E+04 

2.8 

501.9 

5.06E+05 

2.3 

13503.8 

2.19E-04 

4.5 

-1003.8 

4.85E-04 

4.0  - 

■2007.6 

1.00E+13 

0.0 

0.0 

3.00E+13 

0.0 

8604.2 

1.00E+13 

0.0 

5975.1 

1.00E+13  0.0  2629.0 

3.00E+14  0.0  70984.7 
4.00E+12  0.0  0.0 

3.00E+13  0.0  0.0 

3.98E+12  0.0  -239.0 

2.00E+13  0.0  0.0 

2.00E+13  0.0  0.0 

2.40E+13  0.0  0.0 

1.12E+13  0.0  30425.4 
3.98E+13  0.0  61496.2 

1.40E+01  3.1  6931.2 

1.60E+08  1.4  525.8 

1.60E+08  1.4  525.8 

7.11E+08  1.6  7480.9 

9.70E+06  1.9  2361.4 


4.00E+13 

0.0 

4206.5 

5.00E+12 

0.0 

1792.5 

Distribution  A  :  Approved  for  public  release;  distribution  unlimited 
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376.  CH3CH0+0H=CH3C0+H20 

377.  CH3CHO+CH3=CH3CO+CH4 

378.  CH3+CH20H=C2H5+0H 

379.  C2H5+0=CH20+CH3 

380.  C2H5+0=CH3CH0+H 

381.  C2H5+OC2H4+OH 

382.  C2H5+02=C2H4+H02 

383.  C2H5+CH30H=C2H6+CH20H 

384.  C2H6+CH302=C2H5+CH302H 

385.  C2H6+H02=C2H5+H202 

386.  C2H6+02=C2H5+H02 

387.  C2H6+0=C2H5+0H 

388.  C2H6+0H=C2H5+H20 

389.  C2H5+CH30H=C2H6+CH30 

390.  N2H2+0=NH2+N0 

391.  N2H2+0=NNH+0H 

392.  N2H2+0H=NNH+H20 

393.  N2H2+N0=NH2+N20 

394.  N02+0=N03 

395.  NH2+02=NH+H02 

396.  NH2+N0=N20+H2 

397.  NH+0=N+0H 

398.  NH2+0=N0+H2 

399.  CH4(+M)=CH3+H(+M) 

Low  pressure  limit:  0.45000E+18  0. 


1.00E+13  0.0  0.0 

8.50E+10  0.0  6000.0 


1.37E+14 

-0.4 

6596.6 

4.24E+13 

0.0 

0.0 

5.32E+13 

0.0 

0.0 

3.06E+13  0.0  0.0 


2.00E+12  0.0 

4995.2 

3.20E+01 

3.2 

9153.9 

3.00E+1 1 

0.0 

14937.9 

1.70E+13  0.0  20458.9 
3.98E+13  0.0  50908.2 
1.40E+00  4.3  2772.5 

2.20E+07  1.9  1123.3 

1.40E+01  3.1  8938.8 


1.00E+13 

0.0 

0.0 

3.30E+08 

1.5 

497.0 

2.40E+06 

2.0 

-600.0 

4.00E+12 

0.0 

11922.0 

1.32E+13  0.0 

0 

i.O 

1.00E+14 

0.0 

49997.0 

5.00E+13 

0.0 

24481.0 

1.70E+08  1.5 

3368.0 

5.00E+12  0. 

0 

0.0 

2.40E+16 

0.0 

104913.6 

0.90806E+05 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


53 


TROE  centering:  0.10000E+01  0.10000E+03  0.13500E+04  0.78300E+04 


400.  CH3NH+0H=CH2NH+H20 

401.  CH3NH+0=CH2NH+0H 

402.  CH3NH+02=CH2NH+H02 

403.  CH3NH+0=CH30+NH 

404.  CH3NH+OH=CH4+HNO 

405.  CH3NH+02=CH30+HN0 

406.  CH2NH+0=CH20+NH 

407.  CH2NH+0H=CH20+NH2 

408.  CH2NH+0=H2CN+0H 

409.  H2CN+H02=CH2NH+02 

410.  CH2NH+0H=H2CN+H20 

411.  H2CN+02=HCN+H02 

412.  H2CN+NO=HCN+HNO 

413.  CH3+CH3(+M)=C2H6(+M) 

Low  pressure  limit:  0.12600E+42 
TROE  centering:  0.62000E+00  ( 

414.  CH2NH+M=HCN+H2+M 

415.  2H+H2=H2+H2 

416.  NH2+NH=N2H2+H 

417.  N2H2+M=NNH+H+M 

418.  N2H2+H=NNH+H2 

419.  N2H2+NH=NNH+NH2 

420.  N2H2+NH2=NH3+NNH 

421.  2NH2=N2H2+H2 


3.60E+06  2.0  -1192.0 


5.00E+08 

1.5 

-894.0 

1.00E+07 

'  2.0 

6300.0 

6.00E+13 

0.0 

0.0 

6.00E+12 

0.0 

0.0 

6.00E+12 

0.0 

4000.0 

1.70E+06 

2.1 

0.0 

1.80E+05 

2.0 

14800.0 

1.70E+08 

1.5 

4630.0 

1.40E+04 

2.7 

-1609.0 

1.20E+06 

2.0 

-89.0 

2.70E+04  : 

2.0  17300.0 

1.00E+07 

2.0 

4400.0 

3.61E+13 

0.0 

0.0 

I.70000E+01  0.27619E+04 
.73000E+02  0.1 1800E+04 

1.00E+14  0.0  10000.0 

9.20E+16  -0.6  0.0 

1.50E+15  -0.5  0.0 

5.00E+16  0.0  50000.0 
4.80E+08  1.5  1580.0 

2.40E+06  2.0  -1192.0 
1.80E+06  1.9  -1152.0 

5.00E+11  0.0  0.0 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


54 


422.  CH4+H=CH3+H2 

423.  CH4+CH2T=CH3+CH3 

424.  CH3+M=CH2T+H+M 

425.  CH3+H=CH2S+H2 

426.  CH3+H=CH2T+H2 

427.  CH3+CH3=C2H5+H 

428.  CH2T+H=CH+H2 

429.  CH2T+CH2 T=C2H2+H2 

430.  CH2T+CH3=C2H4+H 

431.  CH2T+CH=C2H2+H 

432.  CH2T+C2H6=CH3+C2H5 

433.  CH2S+N2=CH2T+N2 

434.  CH2S+CH4=CH2T+CH4 

435.  CH2S+C2H2=CH2T+C2H2 

436.  CH2S+C2H4=CH2T+C2H4 

437.  CH2S+C2H6=CH2T+C2H6 

438.  CH2S+CH4=2CH3 

439.  CH2S+C2H6=CH3+C2H5 

440.  CH+CH4=C2H4+H 

441.  CH+CH3=C2H3+H 

442.  C2H6=C2H5+H 

443.  C2H6+H=C2H5+H2 

444.  C2H6+CH3=C2H5+CH4 

445.  C2H5+H=C2H4+H2 

446.  C2H5+CH3=C2H4+CH4 


1.30E+04  3.0  8037.0 


4.30E+12 

0.0 

10036.0 

1.00E+16 

0.0 

90607.0 

6.00E+13 

0.0 

15101.0 

9.00E+13 

0.0 

15101.0 

3.01E+13 

0.0 

13512.0 

6.00E+12  0.0  -1788.0 


1.20E+14 

0.0 

795.0 

4.20E+13 

0.0 

0.0 

4.00E+13 

0.0 

0.0 

6.50E+12  0.0 

7910.0 

6.00E+12 

0.0 

0.0 

7.20E+12  0.0  0.0 

4.80E+13  0.0  0.0 

1.38E+13  0.0  0.0 


2.16E+13  0.0  0.0 


4.00E+13 

0.0 

0.0 

1.20E+14  0.0 

6.00E+13 

0.0 

0.0 

3.00E+13 

0.0 

0.0 

2.08E+38  -7.1  106507.0 
1.41E+09  1.5  7402.0 

3.97E+05  2.5  17684.0 

1.25E+14  0.0  8000.0 

4.36E-04  5.0  8300.0 
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447.  C2H5+C2H5=C2H4+C2H6 

448.  C2H4+M=C2H2+H2+M 

449.  C2H4+M=C2H3+H+M 

450.  C2H4+H=C2H3+H2 

451.  C2H4+H=C2H5 

452.  C2H4+C2H4=C2H5+C2H3 

453.  C2H4+CH3=C2H3+CH4 

454.  C2H3+H=C2H2+H2 

455.  C2H3+CH3=C2H2+CH4 

456.  C2H3+C2H6=C2H4+C2H5 

457.  C2H3+C2H3=C2H2+C2H4 

458.  C2H3+CH2T=C2H2+CH3 

459.  C2H3+CH=CH2T+C2H2 

460.  C2H2+C2H2=C2H3+C2H 

461.  C2H2+H=C2H+H2 

462.  C2H2+CH3=C2H+CH4 

463.  C2H2+M=C2H+H+M 

464.  CH+N2=HCN+N 

465.  CN+N=C+N2 

466.  CH2T+N2=HCN+NH 

467.  H2CN+N=N2+CH2T 

468.  H2CN+M=HCN+H+M 

469.  CH2T+N=HCN+H 

470.  CH+N=CN+H 

471.  CH3+N=H2CN+H 


1.40E+12 

0.0 

0.0 

3.00E+17 

0.0 

79349.0 

2.91E+17 

0.0 

96558.0 

5.40E+14  0.0  14902.0 

.05E+14  -0.5 

655.0 

5.00E+14 

0.0 

64700.0 

4.14E+12 

0.0 

11127.0 

1.20E+13  0.0 

0.0 

4.36E-04 

5.0 

8300.0 

1.50E+13 

0.0 

10000.0 

1.08E+13 

0.0 

0.0 

3.00E+13 

0.0 

0.0 

5.00E+13 

0.0 

0.0 

9.64E+12 

0.0 

84448.0 

6.00E+13  0.0  27818.0 

1.81E+11  0.0  17287.0 

4.20E+16  0.0  107000.0 
4.40E+12  0.0  21976.0 
2.40E+13  0.0  -556.0 


1.00E+13 

0.0 

74000.0 

2.00E+13 

0.0 

0.0 

3.00E+14 

0.0 

22000.1 

5.00E+13  0 

.0 

0.0 

.70E+14  -0.1 

0 

.0 

6.10E+14  -0.3  288.0 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


56 


472.  C2H3+N=HCN+CH2T 

473.  CN+HCN=C2N2+H 

474.  NH3+H0N02=H20+NH2N02 

475.  H0N02+0H=H20+N03 

476.  H0N02+H=N03+H2 

477.  H0N02+0=0H+N03 

478.  HNCO+N=NH+NCO 

479.  HNC0+N03=C02+N0+HN0 

480.  HN C O+N 02=C 02+HNN O 

481.  HNCO+M=H+NCO+M 

482.  N20+NC0=C0+N2+N0 

483.  N0+M=N+0+M 

484.  HNCO+CH3=NCO+CH4 

485.  N02+N=N20+0 

486.  CH+N2=NCN+H 

487.  H+NCN=HCN+N 

488.  NCN+N=CN+N2 

489.  C2N2+ONCO+CN 

490.  HNO+N O+N 0=HNN O+N 02 

491.  HNN0+N0=NNH+N02 

492.  HNNO+NO=N2+HONO 

493.  HNN0+M=H+N20+M 

494.  HNNO+M=N2+OH+M 

495.  HNN0+0H=H20+N20 

496.  HNN0+H=H2+N20 


2.00E+13 

0.0 

0.0 

L50E+07 

1.7 

1530.0 

8.00E-01 

3.5  43100.0 

9.00E+10  0.0  0.0 

5.60E+08 

1.5 

16400.0 

1.80E+07 

0.0 

0.0 

2.32E+19 

0.0 

52500.0 

1.00E+12 

0.0  9996.0 

2.51E+12  0.0  26000.0 
1.00E+17  0.0  112000.0 
9.00E+13  0.0  27800.0 

1.40E+15  0.0  148429.0 
1.00E+12  0.0  9935.0 


3.49E+12 

0.0 

-437.0 

2.22E+07 

1.5 

23367.0 

1.89E+14 

0.0 

8425.0 

2.00E+13 

0.0 

0.0 

4.57E+12 

0.0 

8880.0 

1.70E+11  0.0  2100.0 


3.20E+12 

0.0 

270.0 

2.60E+1 1 

0.0 

810.0 

2.20E+15 

0.0 

21600.0 

1.00E+15 

0.0 

25600.0 

2.00E+13 

0.0 

0.0 

2.00E+13  0.0  0.0 
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497.  H0N02+H=N02+H20 

6.00E+13 

0.0 

9800.0 

498.  HON02+H=HONO+OH 

2.00E+13 

0.0 

8000.0 

499.  H0N02+0H(+M)=H20+N03(+M) 

2.47E+08 

0.0  -2860.0 

Low  pressure  limit:  0.68900E+15  0.00000E+00  -0.14400E+04 

N20  Enhanced  by  5.000E+00 

H20  Enhanced  by  9.000E+00 

H0N02  Enhanced  by  5.000E+00 

N03  Enhanced  by  5.000E+00 

NH3  Enhanced  by  5.000E+00 

500.  N03+H202=H0N02+H02  1.00E+12  0.0  8500.0 

501.  H0N0+N03=H0N02+N02  1.00E+12  0.0  6000.0 

502.  N03+N03=N02+N02+02  5.12E+11  0.0  4870.0 

503.  HNN0+N0=N20+HN0  1.00E+12  0.0  0.0 

504.  HNN0+N02=N20+H0N0  1.00E+12  0.0  0.0 

505.  HNN0+N02=NNH+N03  l.OOE+13  0.0  0.0 

506.  HNC0+H0N02=C02+H20+N20  1.00E+05  0.0  0.0 

NOTE:  A  units  in  cm3  mole'cm3  sec"1,  T  in  K  and,  E  in  cal  mole"1. 


Distribution  A  :  Approved  for  public  release;  distribution  unlimited 


58 


